Draft version September 17, 2009 

Preprint typeset using I^'T^]X style cmulatcapj v. 08/22/09 



0^ 
O 
O 
(N 

D 



O 

U 

6 



> 

(N 
m 

0^ 
O 
On 
O 



THE AZTEC/SMA INTERFEROMETRIC IMAGING SURVEY OF SUBMILLIMETER-SELECTED 

HIGH-REDSHIFT GALAXIES 

Joshua D. Younger i-^-^-'', Giovanni G. Fazio\ Jia-Sheng Huang\ Min S. Yun^, Grant W. Wilson"', Matthew L. 
N. AsHBY^^ Mark A. Gurwell\ Alison B. Peck'', Glen R. Petitpas\ David J. Wilner\ David H. Hughes^, Itziar 
Aretxaga", Sungeun Kim*^, Kimberly S. Scott^, Jason Austermann"', Thushara Perera^ James D. Lowenthal'^ 

Draft version September 17, 2009 

ABSTRACT 

We present results from a continuing interferometric survey of high-redshift submillimeter galaxies 
with the Submillimeter Array, including high-resolution (beam size ^ 2 arcsec) imaging of eight 
additional AzTEC 1.1mm selected sources in the COSMOS Field, for which we obtain six reliable 
(peak S/N > 5 or peak S/N > 4 with multiwavelength counterparts within the beam) and two 
moderate significance (peak S/N > 4) detections. When combined with previous detections, this 
yields an unbiased sample of millimeter-selected SMGs with complete interferometric foUowup. With 
this sample in hand, we (1) empirically confirm the radio-submillimeter association, (2) examine 
the submillimeter morphology - including the nature of submillimeter galaxies with multiple radio 
counterparts and constraints on the physical scale of the far infrared - of the sample, and (3) find 
additional evidence for a population of extremely luminous, radio-dim submillimeter galaxies that 
peaks at higher redshift than previous, radio-selected samples. In particular, the presence of such 
a population of high-redshift sources has important consequences for models of galaxy formation - 
which struggle to account for such objects even under liberal assumptions - and dust production 
models given the limited time since the Big Bang. 

Subject headings: cosmology: observations - galaxies: evolution - galaxies: high-redshift - galaxies: 
starburst - galaxies: submillimeter - galaxies: formation 



1. INTRODUCTION 



Though they make up a very small fraction of the local 
infrared (IR) luminosity densi ty, at z <; 1 IR- lum inous 
galaxies (LIRGs and ULIRGs: [Sanders & Mhabel 1996) 
becon ie cosmologically imp ortant (e.g., ILe Floc'h et al.l 
120051: IMagnelli et all 120091 ) and contribut e significantly 



further study is essential to achieve a more thorough un- 
derstanding of the birth and evolution of galaxies. 

A significant population of high-redshift ULIRGs was 
first revealed at 850 /im by the Submillimeter Common 
User Bolometer Array (SCUBA: Holland et al. 1999)_at 
the James Clerk Maxwell Telescope ([Small et al.l 119971 : 



to the diffuse extragalactic IR background fHause r et al 
1998; Kclsall et al. 1998; Arcndt et al. 1998; Dwe k et al 
1998t IFixsen et all 119981: iPei et all 119991 : iDevhn et al 
20091) . They are also thought to dominate the cosmic 
star format i on rate (SFR) dens ity at z ^ 1 ()Blain et al.l 



119991 120021 : IPascale et aIll2009D drive the forma tion of 
luminous qua sars (Sand ers et al.l. 1 988a. b: Hopkin s et al.l 
2006*, '2008b': Coppin et al.' I2008bl : ilvison et al.l 120081 : 
Naravanan et al. 2009b, a) and th e most massive galax- 
ies (Scott ct al. 2002; Blain ct' a ll 120041: ISwinbank et all 
[2006; .Hopkins et al.. .20083; .Viero et alJ 120091 ). As such, 
these ob jects provide power ful constraints on theoretical 
models (|Baugh et ahl 120051 : [S^inbank et al.ll2008[ ). and 



iHughes et allll998l : IsTrger et al.lll998f ). At these wave- 
lengths, the shape of the redshifted far-IR spectral en- 
ergy distribution (SED) counteracts the effect of in- 
creasing luminosity distance to provide an unbiased 
view of du st-obscured star formation out to z 10 
()Blain fc Lo ngair .1993) ■ Since their initial discovery, ob- 
servers have amassed extensive catalogs of t hese submil- 
limet er galaxies (SMGs: for a review, see IBlain et al.l 
[200l) in a number of fields across the sky at both sub- 
millimcter (850/im: Eales et al."1999', '2000'; 'Cowie et aU 
2002; Scott ct al. 2002; Webb et al. 2003; Serjeant et a!] 
20031: IWang et al.ll2004ilCopp in ct al. 200 Q) and millime- 



ter (1100-1200Mm: 
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Ivison et al. 2004 : Greve et all 120041: 
20041 : iLaurent et a l. 2005; S cott et al 
2008; Pcrera etal 



Dannerba uer et aL. . _ 

2006; Bertoldi et al."200 7l:lScott et al 
2008; Wilson et al. 200g nAustermann ct al. 2009) wave- 
lengths. 

However, a more complete understanding of these ob- 
jects has been hampered, in part, by the relatively 
poor resolution of submillimeter cameras (FWHM ^ 
10 — 18 arcsec), which makes identification of multi- 
wavelength counterparts inherently ambiguous. The 
first break-through came with very deep wide-field 
radio continuum surveys, which levera ged the local 
radio/far-IR correlation (for a review, see ICondonlll992 ) 
in combination with statistical arguments to associate 
nearby rad io sources with the submilliimeter emission 
(jlvison et al. 2002. 2007.) . This radio-submillimeter asso- 
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ciation permitted optical spectroscopic foUowup, which 
conf irnied that SMGs he a t high -redshift (median z « 
2.5: iChapman et all 120031 . I2005D and in turn e nabled 
CO s pectroscopic imag i ng (iNeri et al.l[20b 3: Shct h et al.l 
[200l lGreve et al.ll2005l : [Tacconi et aLl[200a i2008'l which 
confirmed that most are young, gas-rich galaxies un- 
dergoing major mergers. While the radio/far-IR cor- 
relation is broadly thought to apply at hi gh-redshift 
(iGarrettI 120021: [Gr uppioni et al.l 1200 3: App leton et all 
12004 iBovle et all 120071 : i Younger et al.ii2009bll . owing to 
the strong dimming of the radio continuum with incr eas- 
ing luminosity distance (see e.g.. lCarilh fc Yunlll999( ) ex- 
isting radio-selected samples are biased towards some- 
what lower redshift 1 ^ z ^ 3 objects. While al- 
ternative counterpart identification t echniques utilizing 
near and mid-IR imaging data e xist (lAshbv et al.l 120061 : 
iPope e t al."2006':'Y un et al.l l2008). these too may be sub- 
ject to biases which are difficult to quantify. 

Since all the above mentioned techniques are inher- 
ently ambiguous, reliable counterpart identification re- 
mains one of the most challenging obstacles to de- 
tailed study of SMGs. Though it currently requires a 
large investment of observing time, this motivates high- 
resolution intcrferometric imaging at the discovery wave- 
length utilizing faciliti es such as the Submillimeter Array 
(SMA: lHo et al.ll2004f) . the Plateau de Bure Interferome- 
ter (PdBI), and the Caltech Millimeter Array (CARMA), 
which provide an order of magnitude improvement in 
the precision of absolute position measurements over sin- 
gle dish instruments and permit unambiguous identi- 
fication of multiwavelength counterparts in higher res- 
olution imaging data. Previous int crferome tric obser- 
vations at millimeter (iDownes et all 1999; Fr aver et al.l 
2p00^'Dannerbauer et al."2002';'Downcs & Solomon"2003l; 
Gcnzel et al. 2003; Kncib ct al. 2005; Grcvc et al. 2005: 
Tacconi ct al. 2006, 2008D and submil limeter ( lono ct alj 
2006i: .Younger et al., ,20°q7l l2008allbl: iWang et al.. ,200/1 : 



Dannerbauer et al.ll2008l : iCowie et al.l 2009f) wavelengths 
have identified unambiguous counterparts for increas- 
ing numbers of radio-detected SMGs, and support the 
radio-su bmillimeter association. In particular, s everal 
groups (lYounger et all ^mft IWang et all IMttI [2?ffl9l : 
iDannerbauer et all 1200^ M. Yun et al. in prepara- 
tion) found that the multiwavelength counterparts of 
radio-dim SMGs provided evidence for a significant pop- 
ulation of SMGs at higher redshift than radio-selected 
samples; a result that has found recent support from 
spectroscopic observations of several individual objects 
jCapak et al.'l2008l; 'Schinn erer et al.ll2008t iCoppin et all 
[2009; Daddi ct al. 2009b, a). 

The existence of large numbers of z ^ 3 ULIRGs has 
profound implications for models of galaxy formation 
and evolution. Recent results from s emi-analytic models 
(jBaugh et al.l 120051 : ISwinbank et al.l [2008,1 have success- 
fully reproduced the observed SMG population at z « 2.5 
- though at the expense of a dramatic departure from a 
universal initial mass function (IMF) in starbursts. How- 
ever, even these tuned models struggle to produce a sig- 
nificant population of higher redshift S MGs (see, e.g.. 
Figure 4 and § 4 in ICoppin et al.l[2009f l. Furthermore, 
given the limited time since the big bang, these high- 
redshift dusty starbursts provide constraints on models 
of dust production in the early universe. 

In order to investigate the nature of these extreme ob- 



j ects, we have exten ded the original unbiased^*^ survey 
([Younger et al.|[2007[ ) of 1.1mm selected SMGs with com- 
plete interferometric foUowup. In this work, we present 
results of high-resolution SMA imaging at 890 /im for 
8 new sources first identified in the AzTEC/COSMOS 
survey ([Scott et al.l[2008) . 

2. OBSERVATIONS 



The COSMOS field ([Scoville et al.l[200l benefits from 
an extraordinary wealth of deep, multi-wavelength cover- 
age from the X-ray to the radio. In this work, we utilize 
i ban d imaging with th e Advanced Camera for Surveys 
f ACS; iFord e"taI|[T99l) on board the EuhUe Space Tele- 
scope to a depth of 27.2 magnitudes (fo r point-sources 
at the 5cr level; iKoekemoer et all I2007D . a variety of 
ground-b ased optical and near-infr ared i maging data (see 
^ani guchi et al.ll2007l : ICapak et a .1120071). ima ging by the 
Infrared Array Camera (IRAC; [Fazio et al.l 12004 ) and 
the Multiband Im aging Photometer for Spitzer (MIPS: 
iRieke et al.ll2004 ) on board the Spitzer Space Telescope 
at 3.6, 4.5, 5.8, 8.0, and 24 /xm to 5a d epths of ^0.9, 
1.7, 11.3, 14.6, and 71 /iJy respectively ('Sa nders et al.l 
l2007f ). and 1.4 GHz radio continuum imaging to a mean 
rms depth of ^ 10.5 /xJy/beam wi th the Very Large Ar- 
ray (VLA: iSchinnerer et al.ll2007n . We also make use 
of photometric redshif ts and stellar mass estimates from 
iMobasher eFall ([2007D . 

The AzTEC/COSMOS survey covers 0.15 deg^ of the 
COSMOS fi eld at 1.1 mm wit h an rms noise level of 1.3 
mJy/beam ([Scott et al.l[200l . The AzTEC/COSMOS 
catalog includes 50 sources with S/N > 3.5f7, of which 
10 sources have S/N > 5(7. For the SMA observa- 
tions we chose the next eight highest significance sources 
(4.5 ^ S/ N <i 5.5) down from the original sample pre- 
sented in lYounger et al.l ([2007| ). This yields, in total, 
an unbiased sample of 15 millimeter selected SMGs with 
complete interferometric followup. Due to the high sig- 
nificance of the sources in our sample, the ex pected false 
detection is <, 0.3 sources (see Fig. 7 in iScott et al.l 
[2008) . 

The SMA observations were performed in the compact 
array configuration (beam size ~ 2") at 345 GHz (full 
bandwidth 4 GHz from the combined sidebands) from 
December 2007 through March 2008. The weather was 
generally excellent (r225GHz ^ 0.08), with typical rms 
noise levels of 1.0-1.5 mJy per track with ~ 6 hours 
of on-source integration. Th e data were calibra ted us- 
ing the MIR software package ([Scoville et al.lll993[ ). mod- 
ified for the SMA. Complex gain calibration was per- 
formed using the calibrator sources J1058-I-015 (~ 3 
Jy, 15° away from targets) and J0854+201 ('- 1 Jy, 
^ 24° away from targets). Passband calibration was 
done using available strong calibrator sources, primar- 
ily 3C273 and Callisto. The absolute flux scale was set 
using observations of Callisto and is estimated to be ac- 
curate to better than 20%. Positions and fluxes of the 
COSMOS sources were derived from the cal ibrated vis- 
ibilities using the miriad software package ([Sault et al.l 

In this context, wc use the term 'unbiased' to refer to a signal- 
to-noise limited sample of millimeter sources, selected entirely on 
the basis of their millimeter emission. While this is not formally 
a flux-limited sample, given some sources with higher fitted fluxes 
but detections below the sigmal-to-noise limit, it is quite close to 
one - especially at the high flux end. 
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Fig. 1. — Multiwavelength stamp images for the new objects iii the sample, including (left to right) the SMA 890/im, Subaru R-band, 
ACS z-band imaging, IRAC 3.6 fim , MIPS 24 ^m , and VLA 20cm data. The red contours over the SMA data represent 3,4, ... X cr in the 
AzTEC 1100 fim map, the red circles over the remaining data indicate the SMA position and are 1 arcsec in diameter - roughly 1/2 the 
SMA beam size, and blue crosses indicate the locations of radio sources within the AzTEC beam. Each stamp is 25 arcsec on a side with 
the exception of the ACS data, which shows a 5 arcsec box indicated by a dotted rectangle. Multicomponent sources are labeled following 
Tabled 



|1995[) . We also incorporated hourly integrations on a test 
uasar J1008+ 063 - which is inclu ded in both the JVAS 
.,Patnaik e t al. 1992VBrowne et al. 1998) and VLB A Cal- 
ibratorjMa ct al. 1998; Bcaslcy ct al. 20Q2,) catalogs of 
compact, flat-spectrum radio sources, and has an abso- 
lute position known to better than 20 mas - to empiri- 
cally verify the phase transfer. 

3. RESULTS 

3.1. Overview 

The multiwavelength data - including high-resolution 
890 fj.m SMA imaging - for the 8 AzTEC 1100/ini- 
selected targets is summarized in Tables [T] and [H and 
Figures [1] and H There are significant ( ^ 3.9a) SMA 
detections for each of the target sources, from which we 
derive absolute positions accurate to « 0.2 — 0.3 arcsec. 
Of those, all but one SMA source are detected in the 
IRAC 3.6 and 4.5 /im imaging (and all but two in the 
IRAC 5.8 and 8.0 fim imaging), while only 5 have MIPS 
24 /im or radio counterparts. We also find (see Figure [3]) 
that those SMA detections with sufhcient I RAC cover- 
age meet the selection criteria proposed by lYun et all 
([2P08). Below, we summarize the data for each target 
source individually. In what follows, we will refer to the 
5/9 SMA sources detected at a peak flux density with 
S/N > 5 as "high-significance" source, while the remain- 
der (4/9) will be referred to as "moderate significance" 



and should be not considered secure detections unless 
they have corroborating multiwavelength counterparts; 
the number of beams in a typical SMA map, in addi- 
tion to experience with coarsely sampled interferometer 
data, leads us to believe that these moderate significance 
detections may be spurious. Furthermore, we make oc- 
casional reference to "power-law" IRAC sources, which 
should be taken to refer to objects whose observed IRAC 
colors are consistent with a significant hot dust contin- 
uum t hought to be powered by an active g alactic nucleus 
fAG N:lLacv et al.|[200llStern et al.|[200alBarmbv et all 
120061: iDonlev et al.ll2007l:lHickox et al 1120071 ). ongoing in 
tense star formatio n (lYun et al. 20081 ). or a complex mix 
of the two ( Youngc ret al.ll2009aD . 

3.2. Notes on Individual Targets 

AzTEC JO95959.34+O23441.O (AzTEC8) - AzTECS 
is detected at high-significance (peak S/N « 10) in the 
SMA image. Its visibility data is best fit by a point- 
source with Fggo^m = 19.7 ± 1.8 mJy offset from the 
AzTEC centroid by 4.9 arcsec. Though the offset may 
seem large compared to, e.g., the signal-to- noise weighted 
error circle (cr » 0.6 x FWHM(S /N)~^ 1.9 arcsec 
for AzTEC8, see llvison et all [20071 ). it is comparable to 
the expectation for an AzTEC s ource with S/N > 4.5 
(4.5 arcsec at 80% confidence, see lScott et al.|[2008[ ) and 
is still well within the 18 arcsec FWHM AzTEC beam 




TABLE 1 

SMA Detections of AzTEC Sources 





Name 


pa 

-^llOOMin 




S/N 


-^890Mm 






AzTEC Offset 






[mjy] 


[mJy] 




[mJy] 


[arcsoc] 


[arcscc] 


[arcsec] 


AzTECS 


AzTEC J095959. 34+023441.0 


5 5+^-3 


19.7 + 1.8 


10.4 


21.6 + 2.3 


0.1 


0.1 


4.9 


AzTEC9 


AzTEC J095957.25+022730.6 




9.0 + 2.2 


4.1 


7.4 + 3.0 


0.4 


0.4 


1.2 


AzTEClO 


AzTEC J095930.76+024033.9 




5.3 + 1.0 


5.3 


4.7 + 1.7 


0.4 


0.4 


0.5 


AzTECll= 


AzTEC J100008. 91+024010.2 




7.4 + 1.9 


8.2 










AzTECll.N= 


AzTEC J100008. 91+024009.6 






10.0 + 2.1 


0.2 


0.2 


2.2 


AzTECll.S'= 


AzTEC J100008. 94+024012.3 








4.4 + 2.1 


0.2 


0.2 


4.7 


AzTEC12 


AzTEC J100035. 29+024353.4 


4 5+^-'-^ 

^■^-1.5 


13.5 + 1.8 


7.5 


12.8 + 2.9 


0.2 


0.3 


1.7 


AzTEC13 


AzTEC J095937.05+023320.0 


4 4+1.3 


8.2 + 1.8 


4.6 


10.0 + 2.8 


0.3 


0.2 


4.5 


AzTEC14.E 


AzTEC JlOOOlO.03+023014.7 




5.0 + 1.0 


5.0 


6.1 + 1.7 


0.3 


0.3 


5.4 


AzTEC14.W 


AzTEC J100009. 63+023018.0 


4 3+1.4 


3.9 + 1.0 


3.9 


4.7 + 1.7 


0.4 


0.3 


6.0 


AzTEClS 


AzTEC J100012. 89+023435.7 




4.4 + 1.0 


4.4 


5.8 + 1.7 


0.3 


0.2 


8.9 



^ Deboosted flux density at 1100 /im llScott et al.|[200a) . The S/N at 1100 was estimated using the peak map signal, rather 
than the deboosted flux density. 

^ The positional uncertainties include both the statistical and systematic uncertainties (see § 13.11 and lYounger et al.ir2007l . for 
details). 

AzTEC 11 shows signiflcant structure in the calibrated visibility data. It is best modeled with two point sources, which we 
designate AzTECll.N and AzTECll.S. In this Table, we separate out the S/N of the signal in the map and the fltted fluxes 
and offsets of the two components. 



(see also overlaid AzTEC contours in Figure [T|) . Of the 
two compact candidate radio counterparts (AzTECS. E 
and AzTECS. W) within the AzTEC beam, only one 
( AzTECS. W; Faocm = 89 ± 11 /^Jy, and coincident with 
a ~ 20 mJy peak in the SMA map) is a significant source 
of submiUimeter emission. AzTECS.E (i^20cm = 139 ±20 
/xjy) is just within the ^ 2a contour but not on a lo- 
cal maximum. Both are power-law IRAC SEDs, and 
AzTECS.E is coincident with a bright MIPS 24 /.tm 
source (F24^m = 820 ± 11)". AzTECS.W is also has 
faintest optical counterpart of the pair. 

AzTEC J095957.25+022730.6 (AzTEC9) - AzTEC9 
is detected at moderate significance (peak S/N w 4) in 
the SMA image. Its visibility data is best fit by a point- 

11 Due to the proximity of AzTEC8.W to this bright 24 ^m 
source ( < 1/2 beam FWHM), it is difficult to obtain a reliable 
estimate of an upper limit on its 24 fim emission. 



source with i^sgoAim — 7.4 ± 3.0 mJy offset by 1.2 arcsec 
from the AzTEC centroid. This fitted fiux density, while 
lower than that inferred from the image plane, is consis- 
tent with its peak flux assuming a point-source structure. 
It is coincident with IRAC 3.6 and 4.5/im, and a radio 
source with F20cm = 59 ± 10 ^Jy, but it is not detected 
in optical, IRAC 5.8/8.0^m, or MIPS 24 ^m imaging 
data. Therefore, despite the moderate significance of the 
detection, because it is coincident with an IRAC/radio 
source we believe the detection is real. 

AzTEC J095930.76+ 024033.9 (AzTEClO) 
AzTEClO is detected at moderate significance (peak 
S/N w 5) in the SMA image. Its visibility data is 
best fit by a point-source with -Fsgo/jm = 5.3 ±1.0 mJy 
offset by 0.5 arcsec from the AzTEC centroid. The 
SMA detection is coincident with an IRAC 3.6-S.O /im 
source which peaks at 5.8/im- though it is statistically 
consistent with a power-law SED - and a faint MIPS 
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TABLE 2 
MuLTi- Wavelength Photometry 



NaiTLG 














4.5^m 






FS, 


^ 20cm 






mag] 


[mag] 


















AzTEC8.E<= 


25 


.69 ±0.11 


25.00 ±0. 


07 


10.8 ±0 


.1 


14.8 ±0.2 


28.3 ±0.9 


69.6 ±2.5 


820 ± 11 


139 ± 20 


AzTECS.W^ 




> 27.8 


> 27.2 




7.4 ±0. 


1 


11.9 ±0.2 


23.8 ±0.9 


34.6 ±2.6 




89± 11 


AzTECQ 




> 27.8 


> 27.2 




2.5 ±0. 


1 


3.2 ±0.2 


< 11.3 


< 14.6 


< 71 


59 ± 10 


AzTEClO 




> 27.8 


> 27.2 




7.1 ±0. 


1 


11.7 ±0.2 


18.5 ± 0.9 


17.3 ±2.3 


114 ± 16 


< 33 


AzTECll^ 




> 27.8 


> 27.2 




29.9 ± 


.2 


42.3 ±0.3 


50.7 ± 1.1 


42.1 ±2.5 


644 ± 11 




AzTECll.N-'' 


23 


.91 ± 0.04 


23.95 ± 0. 


04 














120 ± 26 


AzTECll.S^ 




> 27.8 


> 27.2 
















115 ± 26 


AzTEC12 


25 


.73 ± 0.12 


24.83 ± 0. 


06 


28.5 ± 


.2 


38.2 ±0.3 


59.6 ± 1.1 


56.9 ±2.4 


344 ± 11 


104 ± 14 


AzTEClS 




> 27.8 


> 27.2 




< 0.9 




< 1.7 


< 11.3 


< 14.6 


< 71 


< 33 


AzTEC14.E 




> 27.8 


> 27.2 




< 0.9 




< 1.7 


< 11.3 


< 14.6 


< 71 


< 33 


AzTEC14.W 




> 27.8 


> 27.2 




< 0.9 




< 1.7 


< 11.3 


< 14.6 


< 71 


< 33 


AzTEC15 




> 27.8 


> 27.2 




7.2 ±0. 


1 


12.1 ±0.2 


13.8 ± 0.9 


27.1 ±2.5 


76 ± 11 


< 33 



°' B a nd r' optical imaging data is measured in a 2 arcsec aperture. Upper limits are at the 3-cr level I ITaniguchi et all 
|2003). 

^ IRAC 3.6-8 /x m fluxes are measure d in a 2 arcsec aperture with the appropriate aperture correction. Upper limits are 
at the 5-0- level l lSanders et al.ll2007l ). 

^ MIPS 24 ^m fluxes are measures in a 7 X 7 arcsec square aperture, with the appropriate aperture correction. Upper 
limits are at the 5-(t level l lSanders et al .' 2007). 

VLA fluxes are a Gaussian fit to the imaging data. They do not include corrections for bandwidth smearing, which 
will raise them approximately 15%. Upper limits are at the S-cr level. 

" This source has two candidate radio counterparts with comparable flux density at 20cm: AzTEC8.E and AzTEC8.W. 
However, only one was detected by the SMA. Of the two candidates, one is a bright MIPS 24 /^m source and is not 
associ ated with the submm emission {AzTEC8.E), while the other is not (AzTEC8.W). This source is similar to LH850.02 
from (lYounger et al."2008a'). 

" AzTECll shows signiflcant structure in the calibrated visibilities that are best modeled with a double point source. 
Optical imaging data also suggests a two-component structure, with the northern source (AzTECll. N) having an optical 
counterpart with Zpi-^^t = 2.55. IRAC and MIPS flux measurements are likely a blend of both sources. A double-gaussian 
model has been fit to the VLA imaging data, yielding individual fiuxes for each component. 



source with F24^m = 114 ± 16 /iJy. AzTEClO is not 
detected at optical or radio wavelengths. 

AzTEC J100008. 91+024010.2 (AzTECll) 
AzTECll is detected at high-significance (peak 
S/N « 8) in the SIVIA image. Its visibihty data best fit 
by a double point-source, the two components of which 
(AzTECll.N and AzTECll. S) have Fsgo^m = 10.0 ±2.1 
and 4.4±2.1 mJy and are offset by 2.2 and 4.7 arcsec from 
the AzTEC centroid respectively. There is an extended 
IRAC and IMIPS source in between the two detections 
which is likely a blend of both sources. In addition, 
there is an elongated radio source that is extended in the 
same direction as the submillimeter detection; a double 
gaussian fit to the radio image data yields two sources 
with integrated flux densities of F20cm = 120 ± 26 
and 115 ± 26 ^Jy coincident with AzTECll.N and 
AzTECll. S respectively. Furthermore, AzTECll.N has 
a photometric redshift of Zphot — 1.781 q 23 and a fitted 
stella r mass of log(M*/M0) = 10.9 (jlVIobasher et all 
[2OOI . 

AzTEC J100035. 29+024353.4 (AzTEC12) 
AzTEC12 is detected at high-significance (peak 
S/N « 8) in the SIVIA image. The visibility data is 
best fit by a point-source with i^890/jm = 12.8 ± 2.9 
mJy offset by 1.7 arcsec from the AzTEC centroid. As 
with AzTEClO, the SMA detection is coincident with 
an IRAC 3.6-8.0 /xm source which peaks at 5.8 /xm , 
indicative of 2 ^ z ;^ 3.5 for a starburst (jHuang et al.l 
120041) . in addition to a bright MIPS 24 ^m source 
(^24/jm — 344 ±11 /i Jy) and a radio counterpart with 
F20cm = 104 ± 14 /iJy. 

AzTEC J095937. 05+023320.0 (AzTEClS) 
AzTEC13 is detected at moderate significance (peak 



S/N « 5) in the SMA image. The visibility data is 
best fit by a point-source with Fggo^m — 10.0 ± 2.8 mJy 
offset by 4.5 arcsec from the AzTEC centroid. It is 
not coincident with any optical, IRAC, MIPS, or radio 
sources. There is also a radio source within half an 
AzTEC beam (6.7 arcsec from the AzTEC centroid) 
with an IRAC 3.6 /xm detection and a photometric 
redshift of Zphot = 0.72l°;E5^ (jMobasher et all 120071) 
which is not associated with the submillimeter emission 
- the 890 /zm emission at its location is consistent with 
the noise level of the SMA map. 

AzTEC JlOOOlO.03+023014.7 (AzTEC14.E) - 
AzTEC14.E is a « 5cr peak in the SMA image. 
Its visibility data is best fit by a point-source with 
^890^111 = 6.1 ± 1.7 mJy offset by 5.4 arcsec from the 
AzTEC centroid. It is not coincident with any optical, 
IRAC, MIPS, or radio sources. 

AzTEC J100009. 63+023018.0 (AzTEC14.W) - 
AzTEC14.E is a ~ 4cr peak in the SMA image. 
Its visibility data is best fit by a point-source with 
^890^111 = 4.7 ± 1.7 mJy offset by 6.0 arcsec from the 
AzTEC centroid. It is not coincident with any optical, 
IRAC, MIPS, or radio sources. 

AzTEC J100012.89+023435.7 (AzTEC15) 
AzTEC15 is detected at moderate significance (peak 
S/N sa 4) in the SMA image. Its visibility data is best fit 
by a point-source with i^sgoMm = 5.8 ± 1.7 mJy offset by 
8.9 arcsec from the AzTEC centroid. As with AzTEC8, 
this is wel l beyon d the form al error ra dius estimated by 
both Scott et all 12008) and llvison et al. (2007), and is 
in fact at the edge of the AzTEC beam FWHM. It is, 
however, coincident with a power-law IRAC sources and 
a faint MIPS 24 /xm source (^24^^ = 76 ± 11 /xJy), but 
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Fig. 3. — The mid-IR colors for the sources presented in this 
paper (red stars; hmited to those wi th sufScient IRA C detections) 
as co mpared to those presented in Younger ct al.' (black circles: 
20071) . the radio-selected sample of Chapman ct al. (black triangles: 
20051) . and field sources from the HDFN. The selection criteria 
proposed by [Yun ct al. (2008) arc indicated by solid lines; all of 
the SMGs in the sample arc consistent with these colo r cuts. The 
circle indicates the centroid of foreground sources from lLacv et al.1 
(|200j). Arrows indicate limits. 

is not detected at optical or radio wavelengths. Though 
it only « 4(7 and offset by 1/2 an AzTEC beam from 
the 1100 fim centroid, we believe that detections by 
IRAC and MIPS indicate that the source is real. 

4. DISCUSSION 

The data presented in this work, combined with those 
presented by [Younger ct al. (2007), constitute an unbi- 
ased sample of 15 AzTEC 1100/im-selected sources with 
complete submillimeter interferometric foUowup. Posi- 
tions derived from the interferometric imaging provide 
positions accurate to ^ 0.3 arcsec at a resolution (beam 
size « 2 arcsec FHWM) that is well-matched to multi- 
wavelength imaging data. This, combined with the rich 
multiwavele ngth dataset availab le for the COSMOS field 
(see and IScoville et al.ll2007l for an overview) imme- 
diately enables three investigations: in ij 14. II we present a 
systematic test of the radio-submillimeter association, in 
i) l4.2l we examine the rest frame far-IR morphology of the 
sample, and in ij 14.31 we identify candidate high-redshift 
2^3 sources. In the following a nalyses, we utiUze the 
combined results of this work and lYounger et all ()2007f ) 
for a total sample of 15 1.1mm selected SMGs with com- 
plete interferometric foUowup. As noted in i) l3.1i sources 
at the S/N w 4 — 5 level with no confirming detection 
at other wavelengths may be spurious. This does not, 
however, effect our conclusions. 

4.1. Testing the Radio- Suhmillimeter Association 

The radio-submillimeter association uses 20cm radio 
sources imaged at comparatively higher resolution (20cm 
beam size k, 1 arc sec, absolute astro metric uncertainty 
« 0.1 arcsec: e.g., Ilvison et al.ll2002[ ) to refine absolute 
position measurements for SMGs. This technique lever- 
ages the local far-IR radio relation, combined with sta- 
tistical arguments, to associate SMGs with radio sources 
within the submill imeter beam (see also lPope et al.ll2006l : 
Ilvison et al.l 120071 ) . While this method is efficient and 
physically plausible given the apparent lack of significant 
evolu tion in the far-IR/radio correlation wi th redshift 
(e.g., Ilbar et all 120081 : lYounger et all l2009bD . an inde- 
pendent verification of this technique is important. Our 



Fig. 4. — The expected number of spurious sources associated 
with AzTEC positions for all 15 sources in the sample. We con- 
sider three different definitions of the search area around AzTEC 
detections: (solid line) the full AzTEC beam FWHM of 18 arc- 
sec, (dashed line) the 80% confidoncc separation of 4.5 arcsec as 
estimated by IScott et al. (S08: 2008), and (dot-dashed line) the 
expected positional uncertainty of fitted point-sources to a radio 
map o- f5i 0.6 X FWHM (S/N) ^ (107: Ivison ct al. 2007). The 
number counts were taken from the polynomial fitting function of 
.Bondi et al-i (2008), an d agree with re sults from other deep 1.4 
GHz radio surveys (e.g.. [Hopkins et al.li200 3\ The dotted line in- 
dicates the 3.5(7 limiting flux density of the VLA-COSMOS survey 
(Schinncror et al. 2007, S07:). Even under the most conservative 
assumption of the full AzTEC beam, we would expect at most < 1 
spurious radio source in the full sample of 15 objects. 

unbiased sample of 15 SMGs with direct submillimeter 
positions from interferometric imaging is the first robust 
sample for such a study. 

From a statistical standpoint, we expect the contam- 
ination from spurious radio detections to be minimal. 
Though some observations reveal excess clu stering of 
sub-nijy radio sources on arcminute scales fe.g. jRichardsl 
2000; Georgakakis ct al. 2000), those same observations 
show no evidence of significant anisotropy in 50-100 /iJy 
radio sources on arcsecond scales of the kind considered 
here. Therefore, we can estimate the expected number 
of spurious associations with AzTEC positions assum- 
ing a roughly uniform distribution of foreground radio 
sources and the polynomial fitting form for radio num- 
ber counts in the VLA-COSMOS survey as measured by 
IBondi et all p008). The results of this exercise are sum- 
marized in Figure m We consider three definitions of the 
search radius around AzTEC positions: (1) the full 18 
arcsec AzTEC beam FWHM at 1.1mm, (2) the 4.5 arcsec 
80% confidence interval for high-significance (S/N > 4.5) 
AzTEC detections as est imated via Monte-Carlo sim- 
ulation (jScott et al.l 120081 ). and (3) the signal-to-noise 
weighted error circle of g pa .6 x FWHM (S/N) 
(| Condon! fl997t Ilvison et al.ll2007D . Owing to their rel- 
atively low surface density, even under the most conser- 
vative assumption of a full 18 arcsec search radius, we 
would expect <, 1 spurious association between AzTEC 
and radio detections in the entire sample of 15 objects, 
and under more realistic assumptions none at all. 

A complementary method of estimating contamina- 
tion by foreground radio sources is the P-statistic 
(IDownes et al.l 119861: IDunlop et all 119891 : IScott fc Tou^ 
I1989D . For a given cumulative distribution of number 
counts as a function of flux density nc(> S), one can 
estimate the probability of a random association for a 
source with flux density S and separation 9 as P = 
1 — exp(— 7^06*^). This method has been used by a num- 
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Fig . 5. — The observed submillimeter-to-radio flux density ratio for tlie sample and radio-se lected SMGs (solid triangles: IChapman et aP 
[2005h . The location of GN20, a spectroscopically confirmed z ~ 4 SMG lIDaddi et al.ll2009bl ) is shown as a solid circle. The solid colored 
lines show the range of model tracks from Dale & Hclou (2002), the solid line is Arp 220 (a typical local ULIRG), and the dotted line 
is Mrk 231, which shares many characteristics with Seyfert 1 sources and radio-quiet quasars (e.g.. .Malonev &: Revnoldai200Q ). The left 
panel shows sources in the sample with photometric redshifts (red stars), while the right panel is a zoom-in of the remaining sources 
(dashed lines), and arrows indicate limits. We find that a number of objects in the sample - AzTEC 9, 10, 13, 14. W, 14. S, and 15 - have 
submillimeter-to-radio flux density ratios consistent with a higher median/average redshift than radio-selected samples. 



ber of authors to argue for the hkely association between 
SMGs and radio or mid-I R sources pro ximate to t he sub 
minimeter position (e.g. , 
Ip^ 'Ivis on et al.l[200a 



Hughes et al . 1998; LiU v et al 
Pope et al.l 12006c ilvison et al 



2007). For all the sources in the sample, and a gain as- 
suming the number counts of iBondi et al.l (|2008| ). owing 
to their low surface density on the sky the random as- 
sociation of any significant radio source (S/N ^ 3.5) 
within an individual AzTEC beam FWHM of 18 arcsec 
is highly improbable P < 5%. Therefore, we would ex- 
pect <, 1 total spurious association for the entire sample, 
and far less under more realistic assumptions for AzTEC 
positional uncertainty (see above). Furthermore, consid- 
ering the specific radio detections and their separation 
from the AzTEC centroid, the likelihood of a random 
superposition is even lower, with P << 1%. 

We find that of the 9 sources in the sample with ^ 3ct 
radio sources within the AzTEC beam, there is only one 
instance - AzTEC 13 - in which none of the radio de- 
tections within the AzTEC beam is also detected in the 
high-resolution SMA submillimeter maps; a result con- 
sistent with the statistics estimated above. This target 
has a 8.2 ± 1.8 mJy peak 4.5 arcsec to the north of the 
AzTEC centroid, but no significant detection at the lo- 
cation of the radio source 6.7 arcsec to the south. This 
radio source has a 20cm flux density of 66±10 /iJy, a pho- 



tometric redshift of , 



"phot 



= 0.72 



+0.02 
-0.07' 



and is detected at 



3.6 /im by IRAC. It furthermore has a very high probabil- 
ity of being associated with the AzTEC emission; con- 
sidering its peak flux density and separation from the 
AzTEC centroid we estimate P « 2%. By contrast, the 
submillimeter source is not detected at any other wave- 
length. This is quite surprising, given the frequency of 
bright SMGs with unambiguous position measurements 
and 3 . 6 Mm. c ounterparts (lono ct al. 20061 : 1 Younger et al.l 
[20(r7l I2008at IWang e t al. 2007. 200§!^ Without addi- 
tional information, including deeper radio and IRAC 
imaging, we can only speculate that this source either lies 
at extremely high-redshift and/or has an uncharacteris- 
tically low steUar mass <, IO^Mq (jBorvs et al.ll2005j ). 
At the same time, though we believe this detection is 
reliable, we must also admit the possibility - for the rea- 



sons discussed in § |4] - that an SMA detection at S /N 
« 5 without confirming detections at other wavelengths 
is spurious. 

4.2. The Rest-Frame Far-IR Morphology of SMGs 

In addition to providing accurate absolute position 
measurements, the interferometric imaging constrains 
the submillimeter morphology of bright SMGs. This al- 
lows us to both identify multi-component sources and 
estimate the physical scale of the far-IR. The former 
constrains the contribution of blends to the bright SMG 
population, and the latter is related to the engine driving 
the tremendo us luminosity of these sources (see also the 
discussion m lYounger et'ani2008bf ). 

Our sample of 15 targets contains two reliable targcts^^ 
with multiple radio detections within the AzTEC beam: 
AzTEC5 (AzTEC5.N and AzTEC5.S) and AzTEC8 
(AzTEC8.E and AzTEC8.W). In both cases, SMA imag- 
ing identifies only one of the two radio components as 
the source of the submillimeter emission. In both in- 
stances, the radio pairs have photometric redshifts that 
are marginally consistent with the sources residing at 
the same redshift, and therefore may be physically as- 
sociated; AzTEC5.W (the submillimeter source; see Fig- 
ure 1 in lYounger et all [2007b has Zphot = l-SOto^io and 
AzTEC5.E (« 7 arcsec away) has Zphot = 2.95l?:53 
(more appropriate, these photometric redshifts are not 
inconsistent to within the stated statistical errors); 
AzTEC8.W (the submillimeter source; see Figure [T]) 
has a submillimeter-to-radio flux density ratio consis- 
tent with z « 2 - 3 (see Figure [5]) and AzTEC8.E 
(w 2 arcsec away) has Zp^ot — 2.771q4q. This is con- 
sistent with the statistics of radio sources in the VLA- 
COSMOS survey; under the same assumptions as § 14.11 
the probability of the secondary radio counterparts in 

As noted in § |4] sources S/N < 4 — Scr without confirming 
detections at other wavelengths are possibly spurious. For this rea- 
son, we exclude AzTEC14.W. and E from this discussion, though 
if the source positions are confirmed with more sensitive measure- 
ments this source may represent an example of a truly blended 



8 



AzTECS and AzTECS constituting a random superpo- 
sition of foreground sources are both << 1%. There- 
fore, both objects are similar to LH850.02, a bright 
SMG in the Lockman Hole with two candidate radio 
counterparts of which o nly one was a ssociated with the 
submillimeter emission (jYounger et al.„ .2008a) and sup- 
port the predicted rarity of SMGs arising from confu - 
sion of two lower luminosity sources pvison et al.ll2007D . 
If these two objects are indeed at the same redshift, 
then they necessarily have significantly different far-IR 
SEDs ~ or equivalently, effective dust temperatures - 
with the submillimeter-detected source r epresenting a 
relatively cold starburst (Tdust » 30— 40 K ; iKova cs et alj 
I2006t iCoppin et all l2008al : lYounger et a l. 20G9b), and 
the other a warm starb urst or IR-luminous AGN 
dand ers & Mhabcl 1996; L utz et al.U 998: Risaliti et all 
12000!; [Younger et al. 2009a; Ca sev et al.l l2009). 

Of the remaining sources, all but one are best mod- 
eled as unresolved point-sources, which constrains their 
apparent angular size to ^1.2 arcsec. Furthermore, 
the one source that is resolved by the SMA - AzTECll 
- is best modeled as a double point-source, which is 
consistent with the two components of the extended ra- 
dio emission within the AzTEC beam (see § 13.21 Ta- 
ble [2] and Figure [T|), and its visibility function is in- 
consistent with extended emission such as a Gaussian 
or disk morphology at the 2 — 3cr level. There- 
fore, all the SMA detections in the unbiased sample 
of 15 targets are compact in the SMA imaging data, 
which at z ~ 2 — 3 corresponds to a physical scale 
of ^9 — 8 kpc. While only an upper limit, these 
sizes are consistent with those measured via higher res- 
olution submillimeter (Younger et al.|[20"08bD and radio 
(jChapman et al.l l2004t [Biggs & Ivison '2008') data, and 
rule out cool cirrus d ust (E fst_athiou & Rowan- Robinso3 
120031 : iKaviani et all 120031 ) on > 10 kpc scales as the 
source of the far-IR luminosity in the majority of bright 
SMGs. However, higher resolution imaging will be re- 
quired to measure the size of the starburst region, and 
the data cannot rule out a significant contribution from 
an IR-luminous AGN on sub-kpc scales. 

4.3. Candidate High-Redshift Sources 

Though a number of techni ques have been proposed 
(llvison et al.|[200l iPope et al.ll200l lAshbv et all [2001 
lYun et al.ll20(m " the unambiguous identification of mul- 
tiwavelength counterparts to SMGs represents one of the 
most important challenges to their detailed study. Ow- 
ing to the strong /c-correction at radio wavelengths, this 
problem is partic ularly acute fo r hig h-redshift z > 3 
SMGs. Recently, lYounger et all (|2007) found evidence 
for a population - five of seven AzTEC targets - of 
z ^ 3 SMGs with very faint or no radio or 24 fim 
counterparts. The existence of a large number of these 
high-z SMGs has profound consequences for models 
of ga l axy formation (iBaueh et al.ll2005t iSwinbank et al.l 
120081: iCoppin et all 12000) - which struggle to repro- 
duce these sources, even with the arguably extreme as- 
sumption of a flat initial mass function (I MF) in star- 
bursts - and dust production mod els (e.g.. lGehr3ll989l : 
lMarchenkol200llDunne et al.ll2003D given the limited co- 
ordinate time since the Bi g Bang at these red shifts (see 
also the discussion in § 5 of' Younger et al.ll2007D . Though 
the high-redshift of four individual sources has been con- 
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Fig. 6. — The observed submillimeter-to-24 flux density 
ratio for the sample and radio-selected SMGs (solid triangles: 
IChapman et al.1l2005l) . The location o f GN20, a spectroscopically 
confirmed 2; ~ 4 SMG l|Daddi et al.l [2009b') is shown as a solid 
circle. The grey shaded region shows the range of model tracks 
from|Dalc & Hclou (2002), the soUd Une is Arp 220 (a typical local 
ULIRG), the dotted line is M82 (a local starburst with hotter dust), 
and the dot-dashed line is Mrk 231 (a ULIRG with a signiflcant 
AGN contribution). Red stars indicate sources in the sample with 
photometric redshifts (AzTECS. W and AzTEC12), while the rest 
are shown as dashed black lines, and arrows indicate limits. As with 
Figure [5] a number of objects in the sample have submillimeter- 
to-24 fim flux density ratios indicative of high-redshift (2: > 3). 
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Fig. 7. — Same as Figure[6l but for the observed submillimeter- 
to-3.6 fim flux density ratio and limited to those with 3.6 fim detec- 
tions. Of those with IRAC detections, AzTECO, 10, and 15 have 
flux density ratios indicative of high-redshift {z > 3). However, the 
model disagreement indicates that SMGs have signiflcantly higher 
dust obscuration than is included the model calculations. 



firmed s pectroscopically (ICapak et al.ll2008 | : | Daddi et al.l 
l2009bi lal: ISchinnerer et al.l 120081 : ICoppin et al.l l2009l. the 

role of these objects among bright SMGs, and their 
contribution to the comoving star formation rate den- 
sity (SFRD) in the early universe, has yet to be deter- 
mined. Our expanded, unbiased sample of SMGs with 
interferometric foUowup can provide more powerful con- 
straints, and identifies five additional candidate high- 
resdshift SMGs. 

In Figure [51 we present the submillimeter-to-radio 
flux density ratio for the new objects the sample, along 
with SMGs with spect roscopic redsh ifts and radio de- 
tections from Chapman et al.l (|2005f ). The combined 
effec ts of a negative fc-cor rection in the submillime- 
ter (|Blain fc Longaiil I1993D and a strong positive k- 
correction in th e radio make this quantity a strong func - 
tion of redshift (iCarilh fc Yun|[l999l : I Yun fc Carillill2003 ). 
Those objects in the sample with photometric redshifts 



- AzTECS, 11, and 12 - are all consistent with re- 
sults from radio-selected samples. However, those with 
no radio detection - AzTEClO, 13, 14.E, 14.W, and 
15 - have submillimeter-to-radio ratios that suggest ei- 
ther a higher average/median redshift than radio se- 
lected samples {z ^ 3) or a colder dust temperature. 
However, dust temperatures cold enough to yield high 
flux density ratios (-F890/jm/-F20cm ^ 150) could be ex- 
pected to have extended submillimeter morphologies on 
scales of several arcsec (lEfstat hiou fc Rowan- RobinsonI 
[2003; Kaviani et al. 2003), which is inconsistent with the 
observed sizes of the objects in the sample (see previ- 
ou s section). Th e refore , as with the radio-dim objects 
in lYounger et al.l (|2007f ). we believe these objects are 
compelling candidate members of a population of high- 
redshift sources at z ^3. 

The mid- and near- IR properties of these objects 
are also consistent with a high-redshift. In Figure [S] 
we present the submillimeter-to-24 /zm flux density ra- 
tio, which is also a stron g function of redshift (see also 
iWang et al.l 120071 l2009f ) - though it is important to 
note that the detectability of SMGs at 24 /im depends 
on the relative strength of PAH emission- and silicate 
absorption-features that are redshifted into the MIPS 
24 fim band in sources at z ~ 2, which are quite model- 
dependent. Three of the four high-redshift candidates 
are not detected in deep 24 /im imaging, which is consis- 
tent with z ^ 3. AzTEClO, however, has a ~ 100 /iJy 
24 /im counterpart which, if we assume an Arp 220 tem- 
plate, suggests a redshift of ~ 2 — 3; more in line with 
the radio-selected population. However, this object has 
a 3.6-8.0 /im SED consistent with a power-law, which 
suggests either a significant hot dust contribution from 
either an AGN (i e., Mrk 231) or starburst (i.e., M82; 
see also lYun et al.| [2008). which is consistent with higher 
redshift z ^ 3. The observed submillimeter-to-3.6 fim 
flux density ratio (see Figure [7|), which is also a strong 
function of redshift given a fairly narrow rang e in stellar 
mass and extinction fe.g.. iBorvs et al1l2005D . for these 
five high-redshift candidates is consistent with z ^ 3. A 
full analysis of the near- through far-IR SED, and its 
implications for the starburst properties and redshift of 
these sources will appear in forthcoming work (M. Yun 
et al., in preparation). 

Finally, we can use the observed frequency of high- 
redshift SMGs in the sample - ten total candidate sources 

- to estimate a lower limit on their contribution to cosmic 
star formation at z <; 3. As a crude approximation, we 
assume that the ten high-redshift candidates in the sam- 
ple are uniformly distributed from 3 < z < 5 and have 
a far-IR SED similar to Arp 220, which yields a bolo- 
metric corr ection of LpiR = 2 x (f ggoum /mJy) 10^^ Lq 
(jNeri et all 120031 : lYoung cr et al. 2008b). Given the as- 
sumed cosmology and a Salpctcr (1955) IMF, this yields 
an SMG contribution of SFRDsmg(3 < z < 5) ^ 5 x 
10-3 M0yr-iMpc-3. This is <, 10% of estimates of the 
universal SFRD at similar redshif ts (M adau et al.lll996l : 
iBarger et al][2000t iHopkind l200l . That this is some- 
what higher than previous estimates (|Barger et al.ll2000l : 
llvison et al.ll2002D is not surprising; the sample presented 
in this work contains much brighter objects than typical 
SCUBA 850 /im surveys (median 850 /xm flux density « 
9 — 10 mJy versus w 6 mJy for the Chapman et al. [2005] 
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sample), and sorne have speculated (e .g., iDunlod [20011 : 
llvison et al.) 120021 : lYounger et al.ll2007D that the bright- 
est SMGs lie preferentially at higher redshift. Therefore, 
we find that while bright high-z SMGs are important 
in constraining models of galaxy formation, the high- 
est luminosity millimeter sources (i^noo/^m 4 mJy; 
LpiR. >> 10^"^ ^0) do not dominate star formation at 
3 ^ z 5. 

5. CONCLUSION 

We present results from an extensive campaign to fol- 
low up millimeter-selected SMGs with high-resolution 
(beam size ~ 2 arcsec) interferometric imaging. In 
this work, we have targeted 8 high-significance AzTEC 
1.1mm sources with the SMA at 890/im, resulting in 
six refiable - S/N > 5 ("high-significance") or S/N > 4 
( "moderate significance" ) with multiwavelength counter- 
parts - and two tentative - moderate significance with 
no multiwavelength counterparts - detections. From 
the high-resolution maps, we derived positions accu- 
rate to ~ 0.2 arcsec, in addition to 890 fim fluxes 
and size constraints. When combined with results from 
lYounger et al.l (|2007f ) this constitutes an unbiased sam- 
ple of millimeter-selected SMGs with complete inter- 
ferometric foUowup. From this combined sample, we 
firstly empirically verify the radio-submillimeter asso- 
ciation invoked by previous authors. Second, we find 
that when there are two candidate radio counterparts, 
high-resolution imaging tends to single out one of the 
radio sources as the origin of the submillimeter emis- 
sion, though there is some eviden ce that the two object s 
are physically associated (see also lYounger et aLll2008al) . 
Third, with the exception of one source (AzTECll), all 
of the SMGs in the sample are unresolved by the SMA 
in compact configuration, yielding a maximum angu- 
lar size of <, 1.2 arcsec which, at redshifts typical of 
SMGs, is equivalent to a physical scale of ^ ;^ 8 — 9 kpc. 
Fourth, of the 15 sources in the sample, ten SMGs have 
radio, mid-IR, and near-IR properties consistent with a 
higher average/median redshift than radio-selected sam- 
ples (z ^ 3 — 4 vs. z ^ 2.5 for radio-detected sam- 
ples: see IChapman et al.l I2005D . The existence of such 
a population of high-redshift, hyper luminous starbursts 
has important consequences for models of galaxy forma- 
tion, which struggle to account for such extreme objects 
a z <; 4. 
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